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We report the observation of quantum reflection from a narrow, attractive, potential using bright solitary 
matter-waves formed from a ^^Rb Bose-Einstein condensate. We create narrow potentials using a tightly fo¬ 
cused, red-detuned laser beam, and observe reflection of up to 25% of the atoms, along with the trapping of 
atoms at the position of the beam. We show that the observed reflected fraction is much larger than theoretical 
predictions for a narrow Gaussian potential well; a more detailed model of bright soliton propagation, account¬ 
ing for the generic presence of small subsidiary intensity maxima in the red-detuned beam, suggests that these 
small intensity maxima are the cause of this enhanced reflection. 

PACS numbers: 03.75.Lm, 03.75.-b, 03.75.Kk 


Solitons are non-dispersive and self-localised wave solu¬ 
tions that arise when nonlinear interactions are sufficient to 
overcome the wavepacket dispersion. Since the first obser¬ 
vations in shallow water [1], extensive studies of such soli¬ 
tary wave solutions have been carried out in a diverse range 
of fields, including nonlinear optics and optical fibers [2-4], 
plasma physics [5] and magnetism [6]. In the context of quan¬ 
tum gases, quasi-one-dimensional (ID) Bose-Einstein con¬ 
densates (BECs) may be well described by the homogeneous 
ID Gross-Pitaevskii equation (GPE), a nonlinear Schrodinger 
equation that manifests exact bright soliton solutions for at¬ 
tractive interatomic interactions, taking the form of localised 
density maxima [7]. Experimentally, a quasi-ID limit is ap¬ 
proached by confining the condensate in a highly elongated 
trap with tight radial confinement. While such traps typically 
also feature weak axial confinement, precluding mathemati¬ 
cally exact soliton solutions, the resulting solitary wave solu¬ 
tions retain many characteristics of the ideal soliton [8-10]. 
Previous experimental work has realized both single and mul¬ 
tiple bright solitary matter-waves using ^Li atoms [11-13] and 
^^Rb atoms [14-16], stimulating intense theoretical interest 
(see [10] and references therein). 

Scattering of bright solitary matter waves from narrow re¬ 
pulsive potential barriers has been extensively studied theoret¬ 
ically in the context of matter-wave interferometry [17-22]. 
The nature of the scattering depends crucially on the center- 
of-mass kinetic energy of the solitary wave relative to the 
modulus of its ground state energy. Eor high kinetic energies 
the barrier can act as a beam splitter; the two resulting solitary 
waves can then be recombined to form an interferometer [17], 
with the outcome of the recombination depending strongly on 
their relative phase [23] (as recently experimentally demon¬ 
strated [24]). Eor low kinetic energies, the scattering of the 
matter-wave soliton from the barrier can produce quantum su¬ 
perposition states [25-27]. Previous theoretical studies have 
also addressed the scattering of bright solitary waves from 
narrow attractive potential wells, where the possibility exists 
for the bright solitary wave to undergo quantum refiection. 



FIG. 1. (Color online) (a) Experimental setup. Atoms are cooled in 
a crossed optical dipole trap (not shown), and then transferred into 
an optical waveguide. Additional axial confinement is provided by 
magnetic quadrupole and bias fields. The narrow attractive potential 
is formed using a high numerical aperture (NA) lens to produce a 
light sheet, tightly focussed in the x direction, (b) Absorption images 
of solitary wave propagation in the optical waveguide, (c) Schematic 
showing the position of the narrow attractive potential, relative to the 
trap center and the initial position of the solitary wave. 

Depending on the parameter regime, significant quantum re¬ 
fiection has been predicted for low energy solitons [28], and 
significant resonant trapping when the attractive potential is 
capable of supporting a number of bound states [29]. 

In this Letter, we report the observation of splitting and 
quantum refiection of a bright solitary matter-wave from a 
narrow, attractive potential formed from a tightly focused, red- 
detuned laser beam. We investigate how the fraction of atoms 
refiected varies with the depth of the attractive potential, and 
observe atoms trapped at the position of the well. Surpris¬ 
ingly, we measure much greater refiected fractions than can 
be explained by theoretical predictions for a Gaussian poten¬ 
tial well. We address this discrepancy via extensive theoret¬ 
ical modeling using the GPE, providing strong evidence that 
the presence of small subsidiary diffraction maxima in the red- 
detuned beam, creating a multiple-well structure, is the main 
source of the enhanced refiection. While small subsidiary 
diffraction maxima are generically expected in tightly-focused 
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beams, our experiment is unusual in that they cause qualita¬ 
tive changes in behavior. Our results suggest that carefully 
engineered attractive multi-well potentials may make robust 
beamsplitters for solitary wave interferometry. 

In this work we create stable ^^Rb condensates using the 
method described in [30]. We use the broad Feshbach res¬ 
onance at 155 G between atoms in the F = 2, m/7 = -2 
state [31] to tune the scattering length to positive values, 
avoiding the large negative background scattering length and 
the associated collapse instability [32-34]. Our setup uses a 
levitated crossed optical dipole trap [35] [36], providing in¬ 
dependent control of the trapping frequencies (dominated by 
the optical confinement) and the magnetic bias which is used 
to tune the scattering length (with a sensitivity ~ 40aoG~^ 
around the zero crossing of the Feshbach resonance). We pro¬ 
duce nearly pure condensates of up to 4 x 10"^ atoms at a scat¬ 
tering length of as ~ 200 in an almost spherical trapping 
geometry, (L>x,y,z = x {30(1), 30(1), 42(2)} Hz. The conden¬ 
sate number is reduced to ~ 6000 atoms by further evapora¬ 
tion to facilitate solitary wave production. 

We form a bright solitary wave by loading the EEC into 
an optical waveguide created by an additional dipole trapping 
beam [15] as shown in Fig. 1(a). We first ramp the scat¬ 
tering length to a small positive value, ~ 5 over 50ms. 
We then simultaneously ramp the crossed dipole beams off 
and the waveguide beam on in 250 ms. At the same time, 
we increase the magnetic field gradient to exactly levitate the 
atoms, and ramp the bias field to give a scattering length 
of Gs = -1 Go. This value of Gs minimises the disper¬ 
sion of the condensate as it travels along the waveguide, see 
Fig. 1(b). After the loading is completed the EEC is con¬ 
fined radially in the waveguide beam but is free to propagate 
along the axial direction. The combination of the magnetic 
field gradient, 5', and the magnetic bias field, produces 
a weak harmonic potential in the axial direction, given by 

clJx = 1/2 where ja is the magnetic moment of the 

atoms, and m their mass [37]. This magnetic potential domi¬ 
nates the weak (< 0.1 Hz) optical potential of the waveguide 
in the axial direction yielding overall trapping frequencies of 
Ci>x,y,z = 271 X (1.15(5), 18.2(5), 18.2(5)} Hz. Ey carefully po¬ 
sitioning the magnetic potential minimum with respect to the 
crossed dipole trap we are able to control the motion of the 
atoms in the waveguide. 

We produce a narrow attractive potential well using A = 
852 nm light, focussed to form a light sheet and determine 
the beam waists of = 1.9(2)yum and W3; = 570(40) yum 
via parametric heating of atoms trapped at the focus of the 
beam. The potential well is precisely aligned with respect to 
the waveguide in the vertical direction by mounting the final 
lens in a threaded mount with a pitch of 1.4 jam deg“^. We po¬ 
sition the potential well ~ 22 jam from the minimum of the ax¬ 
ial waveguide potential and release the solitary waves from the 
crossed dipole trap situated ~ 160yum away from the well [as 
shown in Fig. 1(c)], giving an incident velocity of ~ 1 mm s“^. 
At full power we obtain a maximum well depth of 1 yu K x . 
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FIG. 2. (Color online) Solitary wave position as a function of time. 
In the absence of the well (blue triangles) the atoms oscillate in the 
waveguide. With the well present the solitary wave splits, with atoms 
being both transmitted (red circles) and reflected (black squares). 
Lines indicate classical trajectories for free propagation (solid) and 
elastic reflection (dashed). Inset: False colour images taken at (a) 
375 ms and (b) 475 ms with the well present (right) and absent (left). 


In our initial experiment we set the potential well depth to 
its maximum value, release a solitary wave into the waveg¬ 
uide, and track its position by imaging multiple instances of 
the same experimental sequence at different times after re¬ 
lease. Once the solitary wave reaches the well, we observe 
a splitting of the wavepacket and identify three distinct result¬ 
ing fragments: atoms transmitted, reflected, and confined at 
the potential well. We are able to track the center-of-mass po¬ 
sitions of both the transmitted and reflected atomic clouds, as 
shown in Fig. 2. The majority of atoms in the solitary wave 
are transmitted (red circles), following the same trajectory as 
in the freely propagating case (blue triangles), undergoing har¬ 
monic motion in the waveguide (solid line). Around 5 - 10 % 
of the atoms are confined close to the well. The remain¬ 
der of the atoms (~ 25 %) reflect from the narrow potential 
well and propagate in the opposite direction to the transmitted 
component. The turning point of the reflected atoms occurs 
~ 50 ms later than for the transmitted atoms, due to the offset 
of the well position from the trap center. This turning point 
is ~ 20 yum short of the release position, suggesting some en¬ 
ergy is lost during the splitting process. For comparison, the 
trajectory of an elastic collision is shown by the dashed line in 
Fig. 2. 

To explore the effect of the potential well depth relative to 
the energy of the incoming solitary wave we vary the power of 
the 852 nm beam, while keeping all other parameters constant. 
The solitary wave is split and the resulting fragments allowed 
to spatially separate before they are imaged, 475 ms after re¬ 
lease. To calculate the reflection probability, we define three 
fixed regions of the absorption images: transmitted (T), con¬ 
fined (C), and reflected (R), as shown in the inset of Fig. 3(a). 
Taking the sum of the pixel values in each of these regions 
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FIG. 3. (Color online) The percentage (a) reflection (R), (b) trans¬ 
mission (T), and (c) conflnement (C), of atoms as a function of well 
depth for an incident solitary wave with a velocity of 1 mm s~^. These 
percentages are determined using regions deflned in the inset of (a). 
See text for details. 


we define the reflection probability as R/(R + C + T) x 100% 
(values for the transmitted and confined parts are calculated 
similarly) [38]. We find there is no observable reflection from 
the narrow potential well for trap depths <100 nK. Above this 
threshold, the probability of reflection increases sharply [see 
Fig. 3(a)], and the number of atoms transmitted drops cor¬ 
respondingly [Fig. 3 (b)]. For a trap depth of 1 yuK x we 
observe a reflection of ~ 25%. The number of atoms confined 
at the position of the well also increases with increasing well 
depth, as shown in Fig. 3 (c). 

Intriguingly, the observed coefficient of refiection [Fig. 3 (a)] 
is too large to be explained in terms of quantum refiection 
from a Gaussian potential well, 

Vg(x) = -Vocxp{-2x^/{^) , ( 1 ) 

where Vq > 0 and ^ = 1.9 jim. A simple approximate 
argument for this comes from the analytic formula for the 
single-particle refiection coefficient for the similar potential 
y(x) = -Vo/cosh^(x/J) (choosing d ^ £/l.6) [39] 

cos^ 

sinh^(;rM) + cos^ 

where k is the wavevector of the incoming plane wave. Since 
cos^(x) < 1 for all real arguments x, this approximation shows 
that for ^ = 1.9 jum a small incoming velocity (small k) is nec- 
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FIG. 4. (Color online) Reflection coefflcients as a function of poten¬ 
tial well depth for non-interacting wavepackets (a) and bright solitary 
waves (b) in a ID GPE model. Results are shown for a Gaussian po¬ 
tential well with a single potential minimum [yG(-^)]. and for diffrac¬ 
tion pattern [^/(v)] and truncated diffraction pattern [V^^^^^^^v)] po¬ 
tentials with subsidiary minima (see insets). 


essary to observe any reflection, regardless of the well depth 
Vq. For velocities v 1mm s“\ as realized in the experiment, 
this approximation predicts negligible reflection. 

This lack of substantial reflection from the Gaussian poten¬ 
tial Vg(-^) [Efi- (1)] is confirmed by detailed numerical simu¬ 
lations of a quasi-one-dimensional GPE 


dij/(x, t) 


-h- d- 

2m dx^ 


-f y(x) -f U(x, t) + gYo\\l/{x, t)\ 


ll/(x, t), 


(3) 

where U(x, t) represents the time-dependent background po¬ 
tential (see lower curves in Fig. 4). We model the latter as 


U{x, t) = ]-m - xif + - X 2 f], 


(4) 


where xi = -160 pm (x 2 = -22 pm) represents the loca¬ 
tion of the minimum of the dipole (waveguide) potential in 
X, see Fig. 1 (c). The trap frequencies for these potentials 
are ramped linearly over the first r = 250 ms: oJx^O) = 
max{27ryi(T - t)/T,0} and 0Jx2(t) = min{27ry2^/T,2;7ry2}, for 
vi = 30 Hz and y2 = 1.15 Hz. The (static) potential well y(x) 
is centered on x = 0 in these coordinates, and the atoms move 
towards positive x. The nonlinearity gw = 47rNashy_^, where 
we take N = 6000 and y^ = 18.2 Hz. We work with j/r(x, t) 
normalized to unity, and initialize the simulation with j/r(x, t) 
in the ground state of the system for potential U(x,t = 0). 
In agreement with the approximate formula [Eq. (2)], these 
simulations confirm that only very weak reflection (< 4%) 
is expected from the Gaussian potential EgC-^), both for non¬ 
interacting wavepackets [Fig. 4(a), = 0], and for bright 

solitary waves [Fig. 4(b), = -7ao]. We have confirmed that 
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the amount of reflection is not signiflcantly changed by the 
use of a fully three-dimensional GPE model, the inclusion of 
noise in the initial wavepacket, or both of these extensions. 

To qualitatively explain the surprisingly large observed re¬ 
flections we consider the effects of subsidiary diffraction max¬ 
ima in the optical intensity, which occur generically in fo¬ 
cusing optical configurations [40]. Since these are generally 
much less intense than the primary maximum, they are typi¬ 
cally ignored when calculating optical potentials in EEC ex¬ 
periments. However, in the context of our experiment, the 
narrow nature of the subsidiary maxima is potentially sig¬ 
nificant; at least when considered in isolation, they are able 
to produce larger reflection than the primary maximum [see 
Eq. (2)]. Also, the presence of multiple potential wells can 
itself enhance reflection; this is seen, for example, in Bragg 
reflection of BECs from a multiple-well lattice [41] (although 
in our case the rapid variation in well depths precludes a sim¬ 
ilar analysis). 

While the exact structure of the subsidiary diffraction max¬ 
ima in the red-detuned beam is not precisely known in our 
experiment, as a generic model we consider the potential due 
to the intensity pattern of Eraunhofer diffraction from an aper¬ 
ture [40] 


V/(v) = -Vo 




(5) 


and the same potential truncated after the first subsidiary max¬ 
ima; 


T^trunc)/ \ _ j Vj{x) \2xll\ < ^ 2 , 

0 \ 2 xli\>a 2 , ’ 

where a 2 is the second positive zero of the Bessel function 
J\{x). As shown in Eig. 4 (inset) these potentials have a sim¬ 
ilar central minimum to yG(-^), but also one [Vj™^\x)] or 
a decaying series [V/(v)] of subsidiary minima. The results 
of ID GPE simulations for both noninteracting wavepackets 
[Eig. 4(a), = 0] and for bright solitary waves [Eig. 4(b), 

= -Igq] show that the reflection is greatly enhanced for 
both of these potentials compared to yG(-^), the range of 
well depths used in the experiment. The presence of sub¬ 
sidiary diffraction maxima in the beam producing the poten¬ 
tial well thus provides a plausible explanation for the non¬ 
zero reflection probabilities observed in the experiment. The 
similarity of the results for Vj{x) and indicate that 

the oscillatory structure of the reflection coefficient is primar¬ 
ily a transmission resonance effect, attributable to the three- 
well potential composed of the main beam maximum and the 
largest two subsidiary diffraction maxima. 

There are quantitative differences between the experimen¬ 
tal data [Eig. 3(a)] and this generic model; in particular, the 
model exhibits negligible (< 1%) confinement, and oscilla¬ 
tory structure. We have excluded small shot-to-shot changes 
in the incoming soliton velocity due to small (~ +5 jam) shifts 
in the alignment of the experimental potentials as an explana¬ 


tion for the lack of oscillations in the experiment; changing 
the initial displacement of the soliton by up to +5 jam in the 
model leads to a negligible change in reflection coefficient. 
We suspect that these remaining differences may arise from 
effects not captured by our one-dimensional Gross-Pitaevskii 
model, such as the exact structure of the potential well (possi¬ 
bly including time-dependent fluctuations), three-dimensional 
effects, and finite-temperature effects. 

In summary, we have observed quantum reflection of a 
bright solitary matter-wave from a narrow, attractive potential, 
formed by a tightly focused laser beam. Reflection probabil¬ 
ities of up to 25% are measured, with the remaining atoms 
either transmitted or trapped at the position of the potential 
well. Modeling of the system suggests that the exact spatial 
form of the potential well is crucial in determining the amount 
of reflection observed, with the presence of multiple optical 
diffraction maxima, rather than a single Gaussian maximum, 
playing an essential role. These results indicate that carefully 
engineered attractive multi-well potentials could be developed 
as robust beamsplitters for use in solitary wave interferome¬ 
try. In future work we plan to replace the focused laser beam 
with a room-temperature super-polished glass prism (shown 
in Eig. 1), allowing us to explore quantum reflection due to 
the attractive Casmir-Polder potential [42] . 

The data presented in this paper are freely available to 
download [43]. 
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